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The chiral magnetic and chiral separation effects—quantum-anomaly-induced electric cur¬ 
rent and chiral current along an external magnetic field in parity-odd quark-gluon plasma— 
have received intense studies in the community of heavy-ion collision physics. We show that 
analogous effects occur in rotating trapped Fermi gases with Weyl-Zeeman spin-orbit coupling 
where the rotation plays the role of an external magnetic field. These effects can induce a mass 
quadrupole in the atomic cloud along the rotation axis which may be tested in future experi¬ 
ments. Our results suggest that the spin-orbit coupled atomic gases are potential simulators of 
the chiral magnetic and separation effects. 


The recent experimental breakthroughs in generating synthetic spin-orbit coupling (SOC) in both 
bosonic [1] and fermionic gases [2, 3] have opened a new era for cold atomic physics. In these exper¬ 
iments, a pair of Raman lasers induced an equal mixture of Rashba and Dresselhaus SOCs between the 
(psuedo)spin-l/2 internal degree of freedom and the orbital motion of the atoms. Very promisingly, other 
types of SOC, e.g., the Weyl SOC [4-6], could also be realized. The presence of the SOC modifies the dy¬ 
namics on both single-particle and many-body levels and a variety of novel properties have been explored. 
Furthermore, it is very suggestive that the spin-orbit coupled atomic gases may provide ideal platforms 
to simulate intriguing phenomena that have topological origins, e.g., the topological insulators or super¬ 
fluid and Majorana fermions [7-10], the spin Hall effect [11, 12], and the Berezinskii-Kosterlitz-Thouless 
transition [13, 14]. See Refs. [15-17] for reviews. 

In this article, we demonstrate that yet another topological phenomenon, the quantum anomaly, can 
also be realized in a special setup for the spin-orbit coupled atomic gases, namely the trapped rotating 
atomic gases with Weyl-Zeeman SOC. As a consequence of this quantum anomaly, the currents of opposite 
chiralities (see below for definition) are generated in parallel or anti-parallel to the rotation axis. These 
currents mimic the chiral magnetic effect (CME) [18, 19] and chiral separation effect (CSE) [20, 21] that 
are intensively studied in the context of quark-gluon plasma (QGP) produced in heavy-ion collisions, with 
now the rotation playing the role of an external magnetic field. 

The QGP version of the CME and CSE is expressed as 
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for each flavor of light quarks, where ey is the electric charge of quark with flavor /, jy = (ij/yijf) and 
jA - (<A 7 T 5 tA) are electric and chiral currents, N c = 3 is the color degeneracy, and jj and ha are vector and 
chiral chemical potentials. Experimentally, signals consistent with CME and CSE have been observed in 
heavy-ion collisions at Relativistic Heavy Ion Collider (RHIC) [28] and Large Hadron Collider (LHC) [29]. 
In these collisions, extremely strong magnetic fields arise due to the fast motion of the ions [30-32], and 
these magnetic fields induce charge separation and chirality separation in the QGP via CME and CSE 
which in turn lead to special azimuthal distributions of the charged hadrons that are finally measured by the 
detectors. 

It is worth noting that the CME was also discussed in astrophysical context [22] and more recently in 
Weyl and Dirac semimetals [23-27]. In Weyl and Dirac semimetals, the low energy excitations are Weyl 
and Dirac fermions. These materials can exhibit finite /i$ due to their special band structure or by applying 
parallel electric and magnetic fields, and thus open the possibility of realizing the CME. Comparing to 
these previously explored systems, the atomic gases with Weyl-Zeeman SOC enable greater flexibility in 
controlling the parameters and thus provide not only simulators but also a unique mean to exploit new 
features (e.g., those raised by the presence of the Zeeman splitting field) of the CME and CSE. 

Semiclassical approach. — We begin by considering the following single-particle Hamiltonian for spin- 
1/2 atoms (either bosons or fermions) in three dimensions (3D), 

*H = -n- W(p)-cr, (2) 

2m 

where m is the mass, n is the chemical potential, and p is the canonical momentum. The third term ex¬ 
presses a generic SOC where cr is the Pauli matrix and W represents a momentum-dependent magnetic 
field. Its form will be given when necessary. The Hamiltonian (2) represents a two-band system with band 
dispersions e c (p) = p 2 /(2m) — fi- cW(p),c = ±. Correspondingly, we define the chirality of band c to be 
right-handed (left-handed) if c = + (c = -). We note that the chirality we defined here is commonly called 
helicity in literature. We choose the term “chirality” to keep the consistence with the terminologies “chiral 
anomaly”, “chiral magentic effect”, etc. 

Now consider that the atoms are trapped by an external potential V{x) and at the meantime are subject to 
rotation with angular velocity oj which we assume to be a constant. The effect of the trapping and rotation 
can be described byagauge potential - ( Aq,A ) withAo(x) = P(x)-(m/2)(o»xx) 2 -q and A(x) = mcoxx, 
and then the Hamiltonian becomes 

'H = [P ~ o 4(x)] " _ W [p - A(x)] • cr + A 0 (x). (3) 

2m 

Note that the Hamiltonian of atoms with SOC in the rotating frame is generally time dependent and the 
minimal substitution adoptted in going from Eq. (2) to Eq. (3) may not be applicable. However, for situations 
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that, e.g., both the trap and the lasers inducing the SOC are rotating and the SOC is linear in momentum, one 
can find a time-dependent unitary transformation to eliminate the time dependence from the Hamiltonian; 
see Ref. [33] and the references therein. Nevertheless, the present study is restricted to such a situation and 
the use of the time-independent Halmitonian (3) is justified. 

We now derive a set of semiclassical equations of motion (EOM) for the orbital variables x and p. (For 
non-rotating spin-orbit coupled atomic gases, a similar derivation is given in Ref. [34].) At semiclassical 
level, we treat x, p, and cr as classical variables, and their EOM are easily derived from Hamiltonian (3): 


V p £ - V p W • cr, 

(4) 

-V x £ + V X W • cr, 

(5) 

\(T X W, 
n 
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where £(x, p) = [p - A(x)] 2 /(2m) + Ao(x) and W(x, p) = W[p - A(x)]. To proceed, we make an adiabatic 
approximation to the spin dynamics, that is, we treat the orbital degrees of freedom x, p as slow variables 
while the spin cr as fast variable and solve Eq. (6) up to first order in time derivatives of the orbital variables. 
This gives 
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where w = W /W. This procedure is essentially equivalent to solving cr up to first order in h. Inserting 
Eq. (7) to Eqs. (4)-(5) we obtain 


x = V p (£ - cW) + cTi£l px ■ x + cMl pp ■ p, (8) 

P = -V x (£ - cW) - cfiil xx • x - ch£l xp • p, (9) 


where the tensors ft^ e with A, B = x, p are defined as 
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In terms of the kinetic momentum, k = p - A(x), Eqs. (8)-(9) can be recast to more compact and transpar¬ 


ently gauge invariant forms [35, 36], 
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where e c (k) = k 2 /(2m) - fi — cW(k) and 
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are the effective electric field, magnetic field, and Berry curvature, respectively. Equation (14) exhibits 
the equivalence between B and a>. Note that for relativistic massless particles under rotation, the effective 
magnetic field would be momentum-dependent [37], B ~ 2|k|o>, and thus the rotation is no longer equivalent 
to a magnetic field. In that case, the rotation can induce an independent current other than the CME/CSE 
which is called chiral vortical effect (CVE) [38M2]. In the non-relativistic case, the CVE is equivalent to 
CME/CSE. 

From Eq. (11) and Eq. (12) we obtain 

y/Ck-x = V k £ c + eft .E x ft + cft(il ■ V k £ c )B, (16) 

yf(T c k = E + V k £ c x B + cft(E ■ B)ft, (17) 


(13) 

(14) 

(15) 


where \/(j 7 = 1 + c/iB ■ ft, its physical meaning will be clear soon. These are the semiclassical equations 
for the orbital motion of atoms of chirality c. In these equations, the quantum effects are reflected in the 
Berry curvature terms. 

We will hereafter focus on Fermi gases and we will use the natural units ft = k% = 1. 

In the presence of the Berry curvature, the invariant measure of the phase space integration for atoms of 
chirality c needs to be modified to s/G^cPxd'k/dn) 2 [43, 44], With this notification, one is able to write 
down a kinetic equation for the distribution function f c of chirality c. 


d t fc + x • V x / C + k • V k / C = I[f c ], 


(18) 


where x and k are given by Eqs. (16)-(17). In the context of relativistic chiral medium like QGP, similar 
kinetic equation has been derived recently [37, 45—48]. The density and current of chirality c are defined as 


' = Jwfe 
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and the continuity equation, following Eq. (18), reads 
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where we suppose that the collision kernel conserves the particle number for each chirality. If not, there 
would be a term f d 3 k/(2n) 3 \[G^I\f c \ on the right-hand side. The k 0 in the second line specifies the location 
of the Berry monopole which coincides with the band-crossing point determined by W(ko) = 0. The F is 
the total Berry curvature flux associated with the Berry monopole. Its explicit expression is 



This counts the winding number of the map w : 5 2 —> S 2 and thus F e 7i2(S 2 ) = Z. 

If the right-hand side does not vanish, Eq. (21) represents a quantum anomaly for the current of chirality 
c in the form analogous to the chiral anomaly in gauge field theory. This can be seen more clearly if we 
consider a Fermi gas at zero temperature with pure Weyl SOC, W = /Ik, where A is the strength of the 
SOC. In this case, fl = k/(2 k 3 ), V k • fl = 2/rci (3) (k), F - 1, ko = 0, and / c (ko) = 1. Thus the right-hand 
side of Eq. (21) reads c(E • B)/(47r 2 ), which coincides exactly with the 1/(1) chiral anomaly. In this case, 
the conservation of particle numbers of chirality c which is proportional to the volume of its corresponding 
Fermi sphere is violated by the flux of the Berry curvature across the Fermi surface. For a Weyl-Zeeman 
SOC, W(k) = dk + hi with h > 0 being a constant Zeeman field, the Berry curvature is 


n = a 2 


W(k) 

2 IT 3 ’ 


(23) 


the Berry monopole locates at ko = -h/Az, and the winding number F = 1. At zero temperature, we have 
/ f (ko) = 9[pc - h 2 /(2mA 2 )] which vanishes if the Zeeman field h > A Thus the system has quantum 

anomaly only when h < A yj2m/u. The absence of the quantum anomaly when h > A \j2m/j reflects a change 
of the Fermi-surface topology as shown in Fig. 1 . 

We note here that there is no quantum anomaly for Rashba-Dresselhaus SOC W(k) = (Ak y - rjk x , -Ak x + 
T]k y , h) with A, T], h > 0 being constants, because its Berry curvature, 

(A 2 - rj 2 )h 


fi 


2[h 2 + (A 2 + rf)(k 2 x + k 2 ) - 2Ar 1 k x k y ?l 2 


z, 


(24) 


leads to zero winding number. 

Before we proceed, let us comment on the validity regime of the semiclassical approach. To arrive at 
Eqs. (16)-(17) the inter-band transition has been neglected which means that the force acting on the atoms 
cannot be strong. In particular, this require that V|B| <sc VT(p)/d; see Refs. [35-37] for more discussions. 
Obviously, this condition is violated if k is close to the Berry monopole where VF(p) is small. Thus the 
phase space integral in Eq. (21) should be understood to exclude the region |k - kol < A around the Berry 
monopole with A large enough so that we can apply the classical description to particles outside of it. The 
value of A depends on B and E. For example, for pure Weyl SOC, A should be larger than VP; this 
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FIG. 1: (Color online) Upper panels: The dispersion relations e f (kj. = 0 ,k z ). Lower panels: The Fermi-surface 
topologies in (k-, k x ) plane. The green dashed line represents the chemical potential. Blue (red) lines are for right- 
(left-)handed fermions. When h < A y/2m/u, the Berry monopole (the black point) locates inside both the Fermi 
surfaces; At h — A y/2 m/i the two Fermi surfaces touch; when h > A J2mjl, the Berry monopole moves out both the 
Fermi surfaces. 


actually constrains the magnitude of B: because A should not exceed kp, the maximum |B| should not 
exceed k 2 F . This implies that the rotation frequency should not exceed sp in order to guarantee the validity 
of the semiclassical approach. Our numerical simulations will always be within the validity regime of the 
semiclassical approach. 

Chiral magnetic and chiral separation effects. — A consequence of this quantum anomaly is the 
appearance of CME and CSE. To see this, we substitute Eq. (16) into Eq. (20), 
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The first term on the right-hand side is the normal number current, the second term is the (intrinsic) anoma¬ 
lous Flail effect and the last term represents a B-induced current which we denote by j®~ ind : 


j?~ ind = TcB, (26) 

Ac = C J ^^3 (O • V k £ c )/ 0 (27) 

where Xc is the B-induced conductivity (BIC) of chirality c. Let us consider f c to be the Fermi-Dirac 
distribution, 

fc = fc = 8f 1 , for fermions. (28) 

e ps c + I 

In this case, one finds that 

XR = -XL = ^2 ln [1 + ex P (29) 
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FIG. 2: (Color online) The B-induced conductivity, xr • for Fermi gases with Weyl-Zeeman SOC as function of the 
temperature (left panel) and h/(A y/2m/j.) (right panel). Note that xl = ~Xr- 


with eq = n\h 2 /(Imu^ 1 ) - I]. Note that the A and h dependence of^ c is through their ratio h/A, and thus 
when h = 0 the BIC is independent of A as long as it is nonzero. We present the numerical results for xr 
in Fig. 2. It is seen that the BIC is enhanced at finite temperature and suppressed by large Zeeman field or 
small SOC strength. The latter effect is more transparent at zero temperature at which jf = 6{-s c ) and^ c 
can be analytically obtained: 


T c (r = °) = — 2 [ h - 


h 2 

2mA 2 




h 2 

2mA 2 


(30) 


Thus, the B-induced currents disappear when li > A \j2nt/j which again reflects a Fermi-surface topology 
transition (see Fig. 1), in parallel to the absence of the quantum anomaly. 

The above result xr - ~Xl reflects the fact that the right-handed and left-handed atoms have the same 
chemical potential. If the atomic cloud contains domains (a possible realization will be presented in next 
section) in which the right-handed and left-handed chemical potentials differ, say, 


Hr = h+Ha, Hl = h~Ha, 


(31) 


the two BICs, xr and^L, in these domains will also differ in magnitude and Eq. (30) becomes 


c I h 2 

Xc{T = 0) = ^2 U - 


2mA 2 




(32) 


At zero Zeeman field, by inserting Xc(T = 0) into Eq. (26) we obtain the B-induced vector and chiral 
currents: 


•B-ind _ *B-ind , fB-ind 
Jy - Jtf L 


•B-ind _ *B-ind _ |B-ind 
J A ~ JR J L 




(33) 

(34) 


These equations express the CME and the CSE in forms consistent with Eq. (1). 
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FIG. 3: (Color online) A schematic illustration of CME and CSE induced chiral dipole and mass quadrupole in Fermi 
gases with Weyl SOC. First, the CSE drives a chirality separation along the rotation axis with > 0 in the upper tip 
and /.ia < 0 in the lower tip (the left panel). Then the CME in turn drives particle number or equivalently the mass to 
flow away from the center and the atomic cloud acquires a mass quadrupole along the rotation axis (the right panel). 


Finally, we note that the BICs do not receive perturbative corrections from scattering [41,49-52] (See 
the Method.). This originates from the fact that the chiral anomaly is free of renormalization (the Adler- 
Bardeen theorem). 

Chiral dipole and mass quadrupole. — Now we turn to the phenomenology of the CME and CSE. We 
consider a Weyl-Zeeman spin-orbit coupled Fermi gas in normal phase with /./ > 0. Once it is rotating, the 
CSE (34) will drive the chirality to move along B and causes a macroscopic separation of the chirality, see 
Fig. 3 (left). This chiral dipolar distribution naturally forms two separated domains, one with > 0 in the 
upper space while another with ha < 0 in the lower space. Then in these domains the CME (33) in turn 
drives particle number or equivalently the mass to flow away from the center and the atomic cloud acquires 
a mass quadrupole along B, see Fig. 3 (right). In the context of heavy-ion collisions, similar mechanism 
was proposed to generate a charge quadrupole in QGP which may induce a difference between the elliptic 
flows of 7r + and n~ [53, 54] that has been detected at RHIC [55]. 

The above argument is true only at the qualitative level, to reveal the real dynamical process, one needs 
to solve the coupled evolution problem of the right- and left-handed currents. Let us consider a pure Weyl 

SOC and the T = 0 case. The general forms of the right-handed and left-handed currents should contain 

diffusion and normal conducting terms, so they read (the anomalous Hall terms vanish for pure Weyl SOC 
because of the time-reversal symmetry) 

jfl = £p-B - D R Vn R + O-E, (35) 

j l = --J%B - D L Vn L + O-E, (36) 


where D c is the diffusion constant and cr is the normal conductivity, they are linked by the Einstein relations 
cr = DR(dnR/d/jR ) = D/ (d/q /dp/ ). Consider a small fluctuation 5n c in density n c . This will induce a small 
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FIG. 4: (Color online) Mass quadrupole (upper panels) and chiral dipole (lower panels) induced by CME and CSE at 
different rotating frequencies. Length is in unit of 1/ y ma> z . 


departure of the chemical potential from the background value, fi c = + d/C, where ybg = -Aq is linked 

to E by E = Vpbg- Substituting j c to Eq. (21) and keeping linear order terms in Sn c we obtain 

+ TT. T^ B • V6nR - D K y25n R = °’ (37) 

at 4/r z atiR 

f^ B ' V6ni ~ D ^ 5n L = 0. (38) 

at 4n- oiiL 

These two equations represent two wave modes with dispersions Er{ k) = VRk-iDuk 2 , E/(k) = vik-iDjk 2 , 
one propagating along B with velocity vr - B(d/i r / Our) /( An 2 ) and another opposite to B with velocity ry = 
B{dnL/driL)/( 47t 2 ). We call them chiral magnetic waves (CMWs) in accordance with the same wave modes 
found in the context of QGP [56]. It is the CMWs that develops the chiral dipole and the mass quadrupole. 
Unlike the situation in QGP, in our setup the trapping potential will finally balance the driving force due to 
the CMWs and establish a new mechanical equilibrium. The new equilibrium will be characterized by the 
position-dependent chemical potential 

Pc(x) = n - V(x) + A n c (x), (39) 

where V = (m/2) ((lyC 2 + _y 2 ) + oj 2 z 2 j is the effective trapping potential and C 2 = w 2 - a> 2 . The chemical 
potential shift Aq c (x) is determined by the mechanical equilibration condition ,j« = j/ = 0. The numerical 
results are shown in Fig. 4 where we simulate 5 x 10 5 fermions in the harmonic trap. Length is in unit 
of 1 / ymay, and we assume the transverse effective trapping frequency lj ± being kept to be oj- when the 
B-field changes. Experimentally, the chiral dipole may be hard to detect but the mass quadrupole profile 
can be easily detected by, e.g., light absorption images. 

Conclusion.— In summary, we have demonstrated that if the atomic gases with Weyl-Zeeman SOC is 
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1) trapped by an external potential and 2) under rotation, there can appear a quantum anomaly in the chiral 
currents. A consequence of this chiral anomaly is the chiral magnetic effect and chiral separation effect. 
The CME and CSE cause macroscopic separation of chirality and a mass quadrupole along the rotation axis 
in the fermionic atomic cloud which may possibly be detected in cold atomic experiments. 

In the context of QGP, there has been found other transport phenomena that stem from the quantum 
anomaly, e.g., the chiral electric separation effect [57, 58], which may also be realized in spin-orbit coupled 
atomic gases in the similar setup as we discussed in this Report. In addition, it is clear from the derivation 
that the chiral anomaly may exist also in Bose gases. How the quantum-anomaly-induced transports in 
bosonic gases also deserve detailed exploration in future works. 

Method. — Now we show the robustness of B-induced currents against scattering. Let us consider the 
collision kernel I[f c ] in the relaxation time approximation, 

life] = (40) 

T 

where 6f c = / c - / c ° oc B (which is assumed to be small in this calculation) with / c ° the equilibrium 
distribution. We now show that Xc is independent of r. The relaxation time r is assumed to be in¬ 
dependent of the chirality, and we use it to characterize the interaction strength among atoms. Con¬ 
cretely, the relaxation time can be expressed as r = 1 /(nvcr r ) with n the density, v the average velocity, 
cr T = Ana 2 /(I + k 2 a 2 ) the total cross section, and a the scattering length. For fermions, at low temperature, 
r ~ 3nm( I + (kf?a) 2 )/(4k 2 r (k[:a) 2 ). At high temperature, r ~ m 3 ^ 2 Vf fn. 

To proceed, let us turn the E-tield off and assume a steady state with no spatial dependence of f c and B. 
Substituting f c into the kinetic equation (18) in the main text, at linear order in B, we obtain 


Sf c *-—x B) • V k / C °. 

VGc 


The B-field induced current of chirality c is given by 

d 3 k 


•B-ind 

J c 


f ^y[G c : 

J (2n) 3 V 


x(/ c + 8 f c ), 


where the first term gives the result (26) in the main text and the second term 

j3i 


” r / 


' 3 k df° 

—rx(V k e c x B) • V k e c —= 0. 


ds r 


(41) 


(42) 


(43) 
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Thus the B-induced currents are solely given by CME and CSE and are free of perturbative corrections. 
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